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ABSTRACT 


An analytical investigation was made of combined forced 
and free turbulent convection in a circular vertical tube with 
volume heat sources and constant wall heat addition. 

The IBM 1620 digital computers located at The George Wash- 
ington University and also at the United States Naval Academy 
were used to solve the basic equations in which the parameters 
are: Prandtl number Pr; Rayleigh number Ra; friction Reynolds 
number Re*; and the volume heat source parameter F. For fixed 
woes Ol -tilese paramevers -Gne Solution gave the fully developed 
velocity profile, the temperature profile, and the pressure drop. 
The program also solved for the mixed-mean-to-wall temperature 
a@itference, Nusselt number; and Heynolds number. The following 
Veblbaes OL these Paramevers were investigated: Pir = 1, 10,. 100; 
Ra = 0, 34, 54, 104; Re* = 103, 104, 105; and F = .1, .5, 1, 10. 

MIC *GCsultg On Latinar slow problems Were dnvestiegarved Co 
check the validity of the program and they compared favorably 
with known previous results. In the case of pure forced convec- 
tion, laminar flow with a negligible volume heat source, the ve- 
locity profile and the Nusselt number checked exactly with Dr. 


Ojalvo's results for the same problem with no volume heat sources. 
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+ 


The results for turbulent heat transfer in combined forced 
and free convection with volume heat sources showed that an in- 
crease in Prandtl number had a smaller effect than an increas 
in Re* on increasing Nusselt number. The volume heat source 
parameter F had negligible effect on Nu. Hayleighn number had 
no effect on Nu. The pressure-drop parameter C Increased anprox- 
imately an order of magnitude as Re* increased an order of mag- 
mocuge, Prandtl number Pr, Ra, and F had no effect on C for 
Ra < 625, When Ra > 625 and F >1, increasing Pr decreased C¢ 
and increasing Ra increased C; and for Ra + 625 and F <1, Nets 
creasing Pr increased C and increasing Ra decreased C. Decreas~ 
ing F for Ra > 625 decreased C. Increasing Pr had less effect 
than dinereasing Re* on decreasing the mixed-mean-to-wall temper- 
ature difference @,, . Increasing Ra had negligible effect on 
Om . Decreasing F increased @,, . 

The velocity in the center of the tube lowered and became 
less positive as Pr, Re“, or volume heat sources were increased, 
When F ¢ 1, increasing Ra flattened the velocity profile. 

Increasing Re* or Pr flattened the temperature difference 
proiile at we 2erGaclie Ravhacd mecgliiginle effec: on @ . De- 


creasing HW increased (> S00 Mew tare Penoeea tie Jit mei 


ence profile passed through zero. 
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CHAPTER £ 


INTRODUCTION 


me first tncoretical investigations inva combined forced 
and free laminar flow in a vertical tube were conducted by 
Ostroumovlo* and Hallman3, In addition to predieting velocity 
and temperature profiles. Hallman extended the analysis to in- 
clude calculation of Nusselt numbers and pressure drops. He 
considered cases in which volume heat sources were elitner pres- 


ent or not present. He also conducted an experimental investi- 


gation .of the: problem for laminar flow. 
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ied Ceeerevicak JNVesti eeu fale ComoAmcO forced 
eng free turbulent convection in a vertical tube wes conducted 
by Ojalvo and Grosh8, The analysis did not consider cases with 


volume heat sources present. ALI cases of laminar heat trans- 
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fer, pure forced convection or camrbi a Apa free conwmec= 
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ee 
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fu 
a 
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tion, checked exactly with known results. Fr pure forced 
convection turbulent flow were not in agreeiwent with known re- 
sults apparently because too high a value of eddy vigecosity was 
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tSuperseript numbers refer to similarly numbored references in 
fhe Bibliography. 
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results of combined forced and free turbulent convection were 
presented. 

Sackett13 presented data for ccombined forced and free tur- 
bulent convection in a vertical tube with a new relationship for 
the eddy diffusivity of momentum. The results for turbulent 
flow cases indicated some improvement but did not compare Lavor- 
ably with Dr. Ojalvo's results. The most likely cause was attri- 
buted to the reworked computer program. 


Jackson" conducted an investigation to improve tne reiation- 


iy 
S 
Qy 
da 
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Sp for the eddy cilfusivity of momentum. The result 
study are used in the present study. 
The present analysis investigates the transfer of uniform 


thermal energy from the wall of a round tube and uniform ther- 


mal energy from volume heat sources to a fluid ficwing verticaily 


upwards in the tube, Both forced and free convection exist in 
the steady turbulent flow. Only the case of fully developed flow 
is considered and compared with known results, 

im reworking’ the criginal combuver Urcerarm, an crrcer was 
found that would affect Dr. Ojalvo's results for combined forced 
SCs rCe LMP Olont eon vectl@ns Thies ermor dial Mol. occur in 
Sackett's program, thus the apparent cause for the @ifferences 


‘in their results. 
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DiCmanaiy sis nes Dracolcal applicabions in Che fields of 
nuclear reactors and heat exchangers. Its solution wlll show 
the effect of uniform volume heat sources on combined forced 


ena free turbulent convection problems. 
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ANALYSIS 


The problem to be analyzed is combined forced and free 


turbulent convection with untform volume heat scurceés in a 


vertical circular tube whose axis is parallel to the direction 


of the body force. There is to be a net through-flow. 


In add 


Assumptions 


the following assumptions are made: 


lie 


Gr 


Axial symmetry exists for the momentum and heat 
transier., 


All flume properties are constant, except density 
in the expression for body force. A mean density 
Pra” is used for all other density terms. 


Viscous dissipation and axial heat conduction are 
negligible compared with the heat conduction in 
Pie eragla Guwee la cr. 


There is uniform wall heat flux. 


IGLOMevO =UlleeGCoCl D> LOMm-OL auile DEC len 21 ven cabove. 


Tne velocity and temperature profiles are fully de- 


veloped. There are no radial or aneular velocity 
components. 


There is single-phase flor, 


Mie wandre cee aceasta OO ne ee et err a oe en ey Se eS Oe Oe oe BTS 8 Fe ee PH ae eT OEP 2 A err ak & Bn were = Ye 


“Tne Nomenclature is given in Appendix A. 
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(ec VrOW A sGeady, Turbulent and incomprc: 


8, The eddy diffusivity of momentum is given by Jack- 


son's equations. (See Appendix B.) 


Basic Equations 


Pet TL ae he 


ine pasic@equations employed in this analysis 
Continuity, Navier-stokes, and energy equations in 


coordinatesl4, An equation of state is also used. 


On the basis of the above assumpticns and description of 


toe proeplem, the cContinuivcy equation 


oe to NP a Vi) x 
3% (eV) = 0 


reduces to 


ou pa eG) 


dx } 
the Navier-Stokes (Momentum) equations 
dV _& ey + J V7 
= eas D, EF | | aes a fs / 
ea; = Be VP TAVV + ee (Ve) 
reduces to 
Ci oN 2 Val. 20) ( r|. Ms aa) pu) 
Ox UlGe: r 9 Je us f C ar 
2 0 
ends CheMenerey vequarion 
Ceca = Goa. see 
(ree KYO 4 


reduces to 
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(2) 
(3) 
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The equation of state is developed from an expansion of ? in 


a Taylor series in t about the reference temperature t,, and 


Development of Equations 
Equation (1) indicates that 
uszu (vr), (6) 
and equations (3) indicate that 
p =p (x). (7) 
If equations (5), (6), and (7) are utilized, equation (2) can 


be written as 


LAE CEP eke t= Bene © 


Cdr 


The boundary condition cf uniform wall heat flux plus 
the assumotions of constant specific heat, uniform volume heat 


sources, and a fully developed temperature profile require 
at a H (a constant). (9) 
x 


Equation (9) is developed in Appendix C. A new variable a 


is introduced and defined as 


P= Or) St(xr)-t (kK, S) = t-ty - (10) 





In terms of equation (10), equations (4) and (8) become, 


respectively, 

tmepull ~ qh = + Selr[ kre Ce Su] oP) a) 
and 

L dA + Pas GB EO = FPO EES (aay 


Although the density 2 in equation (12) is actually a 
function of x, by considering the selutions of equations (11) 
and (12) to obtain u and @ ata fixed value of x, we can con- 
Beer Pv constant?. our assumption of fully developed flow 


and heat transfer is consistant with these results. Each side 


. 
= 
Ses 
CD 


of equation (12) may be set equal to some constant sin 


have separated the variables-. Thus, 


tas — a ( a constant} (3) 
DFG. 
and. 
| ¢ Emile ; yn g on Sou yO 
rae (PG, + Gaedat) +h Bg 0 =~ Aree an 
C 


The pressure-drop parameter C in these equations takes on the 
value of unity for the special case of pure forced-convecti. 
laminar low, as described by Hallman3. 


10 





In order to solve the set of problems deserilbed, equations 


(11) and (141) are nondimensionalized by using the following 


dimensionless terms: 


CD radius (15) 


We 

d= 16k 04"" Df temperature difference (16) 
eto velocity Gl) 
Ro = Cnn, Oop A DEF 6 pa K Rayleigh number (18) 
= CmnUm Cp Ag" volume heat source (19) 


Equation (19) is discussed in Appendix F. 


Thus, equation (11) becomes 
years J (20) 


% dy vane ae Lm, dy 


and equation (14) becomes 


PROUEIR) = Bf tC 


In accordance with assumption (8), we will use 





: 6 es | 
Es ee Ne ee Sue eae 

e= | <i Ce 
en tre Rai san ‘exp (.01N@ Py) ee 


as given by Lykoudis®. Equation (20) may thus be written as 


A Cle PEI) = roe oe 
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yew - eS rte sed on hme) gmeins 
Boundary Conditions 
ste §46E. BLES Oe a Ponte ORE Ree DOT AL 7S hE Rt ew oe The 


The following boundary conditions are used: 


UG) 8) (24) 


dU = Y’(o) =O (26) 


d@(o%) = 0° (0) aay meer) 


D 


These boundary conditions come from the physical probien, 
Equations (24) and (25) state that the fluid velocity and tem- 
perature at the wall of the tube C7) - 1) are equal, respec- 
tively, to the wall velocity and wall temperature. Equations 
(26) and (27) come from the fact that heat and momentum are 
not transferred across the center line (7 2 Gr Ghewul ne, 


gue sto axial symmetry, resulting im a Zero siope for the tenpe 


ature and velocity profiles at the center line’. 


Discussion cf Equations 
Equation (21), the momentwa equation, and eauation 2a. 
the energy equation, are to be solved by a 1620 IRM digital 


CoMmpuver for U and a as functions oa ». All boundary con- 


ditions are to be satisfied. The value of the pressure-drop 


a2 





parameter C may also be obtained if we use the following form 


Give CONUIiiITy eCqueclone 


Uae ese = — = 


Three dimensionless parameters in our equaticns are the 


> 


(28) 


Prandtl number Pr, the Rayleigh rumber Ra, and the volume heat 
source parameter F. The Prandti number determines the rela- 
tionship between the velocity and temperature dis 
HoOmerilcuspec ta imease vor Fria i S\ boua peoltles are the Gaine. 
hewerci ei number ha =)Gr Pr owiere Une Grastor number ur renre— 


sents a ratio of buoyant forces to viscous forcesl2. Thus, 
free-convection predominates for large Rayleigh numbers and 
for tne special case of Ra = 0, pure forced-convection is reép- 
resented. The significance of variations in the magnitude of 


the volume heat source parameter F is discussed in Appendix F. 


Empirical equations -for < are Piven An A0psne1 2 8, 


ney “are: P a. 
i) = SEE aI een eee PRR 2a A eee 
~ 005 REU-MILHL"G =o025 Re" | por | Res NY < | (29) 
ey eles aa fe es alee 
Sere Tan for 74 SS] a (30) 
éEm— FRe™/F00 . | al 
le 
where - Duw ; {Viction Heynoldads number (32) 
and u* — Mace ew ; Prieta svelociu, (33) 
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Thus, a fourth parameter, the friction Reynolds mumber Re*, 


; : : . WES 
is introduced. In this analysis, when Re* = 0, the ratio 77 


is set equal to zero, thus assisting in checking the validity 


of the program by comparison with laminar flow problem results. 


Beg el ee ee) 
The solution of equations (21), (23), and (28) will yield 
U and © as functions of n for a given set of values of the 
parameters: Ra, Pr, Re*, and F. The pressure-drop parameter 


C Will also. be obtained. 


These results may be extended by calculating the rollowing 


useful quantities: the dimensionless mixed-mean-to-wall tem- 
perature difference @,, ; the Nusselt number Nu; and the Rey- 
nolds number Re, 

Haliman3 gives equations for the first two of these quan- 


Gites... theycare. 
| 
10) es ee [ ONO iig (34) 


and 


Na Bie Fn Pa (35) 


see Appendix D. 


Anvexpression for the Reynolds number as derived in 
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RopenGis i. 


ee = 


+ 


— 


The result is 
Re” )S 
8C + Ra Dm Zt F 





(36) 


Where Une = Siem 1eetor upward alow ane tne = sign ae for 


net downward flow, 


2 


— 





CHAPTER IIL 


METHOD OF SOLUTION 


The method of solution is based on the use of a digital. 
eompuver. This method will decrease congidgerably une amouns 
of hand calculations required in the solution of our problems. 
The possibility of errors are thus reduced. This method of 
solution will also allow us to investigate many ditferent 
problems and compare some of these results with known analy- 
tical and experimental data. 

The independent parameters used are Ra, Ne. anand “i. 
Hie havle teiununMbel Gs “HOSS GO Measure Tie exUenu. OF Une 
free-convective effect on the flow. The friction Reynolds 
number determines the value of “_ as a function of 1 : 

The Prandtl number determines the value of 6 . The volume 
heat source parameter is chosen to deseribe the thernal ener- 
gy convected downstream to the heat generated in the fluid. 
With these input quantities fixed, the method of solution pro- 
eeeds as follows: 

Le Ue oe cl a a (37) 
TS assumed ose sini Toless FOr Ge Velocity morte.  Ahis 
equation satisfies the continuity equation Lundy --~ ie 


16 





and the boundary conditions U(1) = 0, and U'(@) = 0. It is 


the limiting case of the parabolic profile in pure forced- 


convection laminar flow. 
©, Equation (23) is integrated with the aid of boundary 
= oa to obtain | 


war 


condition or o 
a: | EER Jue 7) ak ( 
= cd — : 38) 
Q NII beFr sy [] to" Pr Sx 


Equation (38) is integrated using boundary condition 


Sr 
@ (1) = O to obtain 
(39) 


p= fin 


as a function of 1 
Equation (21) is integrated using boundary condition 


y i 
ut (2) = 0 to obtain 
/ Re, : 4 ©) d as . 
Be a et et ts 
En , 
ies 8, | [ u 
The pressure-dropv parameter C must have an outrut numer 
f pure 


deal value for the computer. For the limiting case o! 
forced=-convection laminar flow ¢ = 1. 
S integrated subject to boundary condi- 


5. Equation (40) 4 


tion U(1) = 0 to obtain 
4 
=f" "Udn (41) 
iS romiact on of mal j 
6. The continuity equation (28) is checked to see if 
equation (41) satisfies it. If Ue does not satisfy this equation 





the value of C is changed until it is satiseried. 
7. The function Up is checked against Uy by comparing 
U'o(1) with Uap iGr) to see if they are within 0.1% of each other. 


8, Equations 


z( Re’) 


ene, Re@= Fiera) 


on 9 Re Ora A FF 
are used to calculate O;,, Nu and Re, respectively, if the 
check in Step 7 is satisfactory. The results are then printed 
out as given in Step 10 below, 
9. If the check in Step 7 is not satisfactory, Uy is used 
dome ibitial assumovilon,- repiaciwns UW anm Coins through tne 
PEOCCOULTE Aeain, Starline at ~oecD.e. 


107. Pring out 


aoe 
Problem Number } 
Ra, Re Pe Cr 
24 Op Uy 
: ! , 4 
a O O O 
C De Nu nile 
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Uy and U; are both printed out to ensure that the criterion 
POr -CWeCkKkine Chem 1 DueD / 25 Valid, 
These steps may be shown in a simplified flow diagram for 


wie =calculavion, Fagure 1, 


a 





CHAPTER IV 


COMPARISON OF RESULTS 


Table 1 lists the following calculated results: 0, 04, 


Re, and Nu for all problems solved by the 1620 digital computer. 


TABLE, 1 


RESULTS FOR ALG COMPLETED PROBLEMS 


Input Parameters Calculated Results 
Problem 
Nol tia) oie ~ ies ny C Dr Ne Nu 

he baa Ol ll CumOl COO CONY CC OLN 
kO3B OO Of 2220 eo 1.0 -915.85 ie 36 

Be eheaniiel fT Vows WeOnoinicd Forced oma efiece 
hoor 8&1 Ol e225 70 ESS Ee == 3.74 

Can Hirbulent tio. pure Lorceduconvceul on 
ao 0 103 1 5 Ta 3 ,036 17000 55.4 
380 0 102 1 sO Tee OOo! 17000 0.0 
383 OF alee iy Ono) (Aes: sar 17000 50.3 
385 Or AOS 10 5 We8 youl 17000 196.0 
peo 3G 10>’) 10. 170 9 47.3° —=.00052- 417000 0.0 
389 Or sle> Ou O40 ie ou 17000 185.0 
391 Opole | 100 5 7.3 ,002 17000 983.0 
302 On 10-4) dloon “ec 78 - 000064 17000 0.0 
395 Opa 10> |) TOCn son Wee ~ 038 17000 952.0 
397 0 10% i 5 53.9 E005 232000 400.0 
398 Oo 10% is SIeOM Gone -,00017 232000 0.0 
l00B O 300 iY ELOEO 230° = 207 1030 17 
Lol Oo 104 le) ho eOne NS oo 005 232000 Asis (6 
103 0 104 10 Sr) See9 .0011 232000 1750.0 
toh o 104 10 1,0 53.9 +,000019 232000 ome 
ogee OMe ® IO .L1O.Om Goro) 405 232000 1720.0 
Log ©. Os +) 106 5 See .00025 232000 7900.0 


ND 
© 





Table 1 (Continued) 


Problem 
No’. 

Gs 
H1O O 
413 O 
AVS O 
LLG O 
119 O 
sei, 0 
ee O 
Los O 
4. O7 O 
1.28 O 
429B O 
55 81 
56 81 
57 Sl 
60 81. 
61 81 
62 B81 
63 81 
66 out 
67 81 
68 81 
69 81 
72 81 
73 B31 
74 81 
ie 81 
78 81 
79 81 
80 81 
81 81 
8} 1 
85 81 
86 Sl 
87 81 
90 Sl 


ha 


Input Parameters 


Calculated Results 


Nu 


OF 
NSO 
S502 

Oz 
3040. 

11700. 

0. 

LUE LO(e). 
18600. 
or 
386, 


Doe 
DoE 
QO, 


5O; 


Re® Ppp F C Om Re 
Turpulent flow, pure forced convection (Continued) 
LO 100 1.0 53.9 -.000002 232000 
Hoey SCO MONO. Sano 0046 23 2000 
o> 1 5 430.0 .00063 2900000 
1.02 1 1.0 430.0 -,000021 2900000 
102 i 10 Or se OrO One 2500000 
10 10 55 430.0 .00017 2900000 
102 10 1.0 30.0 ~-.0000034 2900000 
102 Tom E1OSOn *)00n 6) = .00nn EGOC000 
102 +100 5 430.0 0001 2900000 
10? 100 1,0 430.0 -,0000014 2900000 
10+ ios 53.9 10.34 23.2000 
De ter boulent low. eonltmocen Orced amd alree 
103 1 pall 5.8 .067 16600 
103 1 us eve £036 17000 
tO eee aca OOLG 17000 
102 EO Oen NTeas nana 17000 
103 10 il 6.9 .018 16900 
162 Sto : 5 a on 17000 
162 Tope 7.3. ~.00054 17000 
102 10 10.0 ee iol 17000 
102 | 106 = 72 0057 17000 
02) Bator Sieben Coca. Oe 17000 
Oe PCO alan yeycs eee oVorclersyl 17000 
10; HO OHNO. Onn e=Oai: 17000 
10, 1 16 53.7. ~»=©.0091 232000 
10, H 5 53.9 005 232009 
10) i te 5e-6" -=Goets 232000 
10 I 1OnO) 5360) 27058 232000 
10) 10 mi Sere. ,002 232000 
10), 10 aS Soe? OO 23 2000 
10) 10 1.0 53.9 -,00002 232000 
10), 100) 1ONOm 52-50) = .051 232000 
10, 100 1 53.9 =, 00045 232000 
10, 100 5 53.9 .00025 232000 
10, 100 1.0 53.9 -,.000002 232000 
Of) woo GG: oe Ane OOl, 23 2000 


fas 


@ VON s) (Voor e) Gone) CVG Ss eles 6) ore) one: 


ee) 2 ONO ae 


tj 





Table 1 (Continued) 


Problem 
No. Ra Re* Pr I ¢. 
D. Turbulent flow, combined 
Ql 81 102 i ieee erate 
92 81 10? 1 5 430.0 
Gao) e102 i ae Sect 
96 81 10> iO O me Ono 
97 81 102 #10 .1 427.0 
98 81 102 10 °&.5 430.0 
99 8), ae? 1G io 186.0 
102 81 102 10 10,0 430.0 
103 Si eae> 100 Si shea? Ne 
104 81 102 100 5 430.0 
105 81 102 100 1.0 430.0 
108 81 102 100 10.0 430.6 
164 625 102 1 Bee GO 
165 625 10° (eeeieOn bak 
169 625 102 10 me Peo 
170 625 102 10 6 I GO 
171 625 102 TOm, sO (me 
174 625 102 i weatne eo 7.6 
M5) 66250102 100 1. 6.7 
fim6 | G25 lo 100 65 7.2 
iy 9625102 100 1,0 7.3 
180 625 103 100 10.0 7A 
181 625 10; it Ald PUSahe 
182 625 LUO 1 oS. Sen 7 
183 625 Os i INO Csero 
187 625 10° 10 wy, bSee 
Mes. Noe5mi0° (10° 65. 63.8 
MQ Noes m0! “To. 10) Sacwo 
192 625 Oe TO LOMO Mince o 
MmO3, "625-9110" 100 1 53.8 
Mol 625 Blot #1000095). 53.9 
15 UMCe> EO!  1OOe 1.0 5a Oo 
198 625 10! 100 10.0 53.9 
199 625 10° 1 .1 427.0 
200 625 102 1 PS BOs 
201 625 102 ip. stioe eae 
205 625 102 10 .1) 427.0 


Input Parameters 


fe 


forced and 


Calculaced Results 


Om 


OOn: 
00063 
OOeCeL 
SOL 2 
,0003 
POOo1/ 
.00000345 
,0031 
.O0019 
,O001 
.000001 4 
.OO19 
036 
J0035 
JOLO 
FO1O 
~00056 
194 
,003/ 
O02 
,OC00S4 
035 
{0091 
005 
.00017 
OO? 
,OO1 1 
.000019 
AGIA 
,OOOK6 
,00025 
,000002 
0046 
OO: 
.00063 
,000021 
0003 


Re 


Nu 


free (Continued) 


290000 
24500000 
2900000 
2900000 
29200C0 
C3IOO0O000 
CGQOO0O 
E909000 
2920000 
29Y0GOO0O 
2900000 
2900000 
16700 
1 (O00 
16200 
17300 
1 /OOO 
32000 
16900 
VAL Oe 
£7000 
18800 
PEALOLON®) 
232000 
232000 
232000 
232000 
232000 
23 3000 
232000 
232000 
232000 
232000 
2920000 
2900000 
2900000 
2920000 


S120 20 

Se OR 
0.0 

3040, 


YO 
(ee) 
oe, 
OO OPO OOO 1 OOOO 6 2 © O7Ore.© © 10 OO 2'O1o © G2) 





Table 1 (Continued) 
Input Parameters 


Calculated Results 


Problem 
No. Ra Re* Pr F C Om Re Nu 
i a. oo forced and free (Continued) 
. 5 = 5 5ay Oro FOCOI 2OLOOOO} 117002 © 
ao. i _ 1.0 ge ~-,0000034 2910000 O20 
mec. o. G.. WOO SO 10 see OOs) 2910000 11500.0 
le erate) 0001 
- p 9 2920000 18600.0 
5 10 100 15. 2180n0 0001 6 
Meee is tec 2PILOOOO 5 +416 50020 
ee ioe. Gee a Coca ~,0000014 2910000 O20 
a fea 103 : e i202) =. O01 2910000 18300.0 
aH 1 303 : * 8.3 ~.0035 17200 O26 
a Pe a3 = .O 7 eS) 27600 see Ti 
- at ee ™ ae 1. now 15700 196.0 
=e ea! abe - = (ES k= 20005)! 17100 0.0 
-. a a EO. “Ueno 3204 ley 4@le) 185.0 
od aoc ee a Cee 002 16700 984 .O 
ie 1 Fe a fae t= OOOO72 17000 O20 
18 rol yolk . Cet OSs 17300 948.0 
aie oe pe : ; Ot 0091 223000 595.0 
i OH rod : ae Sub el 005 230000 400.0 
ae oe as : ue =e -,0002 231000 Cro 
aS rot rol 7 : S ie O04 232000 383.0 
= rol roll me ol Wen 002 230000 170.0 
a sod ol an so “DS ae OO11 232000 1750.0 
C4 on - 1.0 53.9 --.000019 232000 0.0 
. oe me aa Sao: 702i 232000 1720.0 
$8: rod ot : 52.5 OO0O45 231000 7890.0 
a so! ro! oe Sy oes .00025 232000 7910.0 
i ie pe sh 0: “5249. = 000002 232000 0.0 
10) a. 4 TO sho 54 Or = OO1G 23 2000 7850.0 
a a : fa ine aout 2920000 3120.0 
os a as : se earn V0063 2900000 3150.0 
i : x0 130.0 --, 00005 2900000 O50 
Pons ale i” O20) FLSO%G, 207 
367 10% 105 10 1 4 : A be tine 
a =p: 05 “1 26.0 00003 2920000 11700.0 
10 Pare 0 40 melelent f 
mo oe me : 7 2910000 #£11700.0 
i me et . Pe pen -,00000344 2910000 On0 
is + Sige - pa —iOOst: 2910000 £11500.0 
a iy os : 127.0 .00019 2920000 18600.0 
no 100 25 430.0 0001 
ee a ; 2910000 #£18600.0 
ah o oc en) 130.0 =sOOOOCO1H =2011C000 0.0 
| 100, CUCrO. Wao a2 Oene 2910000 18400 .0 





Part A consists of the problem in laminar flow, pure forced 
convection with negligible volume heat sources. Fart B con- 
Sous Of a provlem in laminar 110n,” Combined 1To0orced ard Tree 
convection with volume heat sources. Part C consists of 28 
PeOolens JiwelMrbubecne Oval pure TPorced <onvecbhion with Volumes 
heat sources. Fart D conSists of 100 problems in turbulent 
MWCO UoLiiet 1 OPCee ane Tree GONVeECLLOl Wit volume Neat 
SOUrCCS« 

The IBM solutions to all problems are an attachment to 
tae presenc analysis, Copies of tne fortran program, tne data 
Cares.canad the SoOlucions Ge all problems on 18 cards are 1 


mie phic s  O: Dige Ojalve 


Denes Lom 


aera 6 2d se Ss ees 





The velocity profile and the Nusselt number of the laminar 
flow problem in part A, Table 1, checked exactly with the re- 
sults of Ojalvo® and Haliman>, A negligible volume heat source, 
PoLOOO; Woaseticed Am Eials Reopen ©O Checle Tue -valiaity of “the 
program. The leminar fisw proodlem, combined forced and free 
CONVCCGlOU Wo Youle Ties. Sources @f. part 2s Table Liatecon. 
pared favorably with Haliman. It also compared favorably with 


Ojalvo even though no volume heat sources were present in his 


analysis. 
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Problem 397 was chosen as typical of the problems in part 
C. This problem was compared with test problem /29B and the 
universal velocity and temperature profiles of Eckert and 


Drake“. 


The values of U and Oren problem 429B with negligible 


wolume heat sources, FilO00, Were checked by using the defini- 


; 2 
tions 
ut 


Hl 


UTI, (112) 
+- 
Y 


cat)" = inset ae ee 


7 a dy) 
as Vw Pm Cy eee ™ 


with the definitions of the various dimensionless quantities 


[ 


wae ) (113 ) 


to derive the comparison equations 





uw am (U) / (45) 
Re 
yoo = phen , (46) 
(att)- 1 PrRe™ (6) | : 
and © Ge Oa ( ) (17) 


Figure 2 shows the comparison of ut and SEs) for prob-- 
lem 429B with the universal velocity and temperature profiles. 
The values of velocity coincide almost exactly up to yt =100, 
then they start to diverge slowly below the universal profile. 


The values of temverature fall slightly above the universal 





- ~ +: sy 994 ‘ bs ot eek, ae ee Beene . ane Ye 2 naa : 
tempe rature profile, 4 QS Velues ofF problem sol 7 eoncide exac pay 
. ~ ae 


é 


ie ak Ao" a lee ay ee. nie eat. Jobe ry ee ees Pe eee Le C 
with those of problien 4295 excaent tne temperature starts to 


as ee 
'=300., Volume 


Poa 
bt 
ee 
zy 
A ted 
MN 
a: 
Eri 
Q) 
cr 
ree 


diverge above the values 


@ 


Mean sources lad Mociiei pia crilGoe on tne Velceiluy and wenper— 


ature prortles. 

The differences in tne values of these two problems and 
Peeples Ore Oly ere ero iyem oO Che Moro ved express 
sion for the relative viscosity. A peered. eee 2 COW OUOn Ge 
Oe that is an increase in the volume heat sources, had very 


Littie efreet on the vrofiles. 

All of the problems in part C with a Prandtl number of 1 
aM — and With a voiume heat Source paramever F=lO are plotted 
in Figure 3, which shows the varlhation of Nusselt number with 
ian Hower. Whee OP es ea ores equation! Lom sure Lorcea- 
convection turbulent heav transfer is alse plottad, In both 


a 


Prandtl number comparisons, the results of the present analysis 


é 


Comvare almuis Gxactly with tne curcves of Ghe Diccus-Psei ter 

equation. Volume heal soucgsvon wad mee eile efrecy on the 
HiaCLOnvOL NUGSSL Tuner WL bh Neville mienber, 

Pitures 4 through 10 shoy resnlss for tae problems in tux 


bulent flow, combined foeus! and Free convection. The Nusselt 


Lan 4 ieee i eee TORE AWS, tn eae Ce CCT oe ay 
PoC r? 21S DLGCted CO ia eM Cr at CNN Ce Ce es ik qihe 


5 





Nusselt number is increased less than an order of magnitude as 
Re* is increased an order of magnitude. Increasing the Prandtl 
number has less effect on increasing the Nussult number than 
increasing Re*. Increasing the Rayleigh number had no effect 
She WusselLt number, -Vecreasinge i nada — tVEtle Vicressings 
effect on Nu. 

several test problems were run to determine the limitations 
CEAcLNe COMpUucter Lor Re®=102, No solutions were obtained on the 
IBM 1620 computer for Ret=10-, P=.,1 or less, and Raz625 or more. 
Only by decreasing Ra or increasing either Re* or F would the 
computer give a solution. A large Ra in combination with a 
small F (a large volume heat source) would physically appear 
TO be opposing each other. Convective forces predominate in 
a large Ra, whereas with a small F, considerably more heat is 
being generated than is being convected downstream. 

Figure 5 shows the pressure-drop parameter C plotted 
agalnst Rayleigh number. C increased approximately an order of 
magnitude as Re* increased an Beeler Of Macmbtuae., for hayveiah 
number less than 625; Ra, F, and Pr had no effect on Cc. For 
values of Ra over 625 and F?1, increasing Pr decreases C and 
increasing Ra increases C. If F<l1, increasing Pr increases C 


ands nereasine na Wecredses Coe iese (curves. Slow Ula Goeianeces 





slope for a value of F approximately equal to 1. For a given 
he’ eane Pr “he point av Which C changes slope is a function of 
Ra. For small values of F and Re*, and large values of Ra, C 


would appear to go negative verifying Hallman’ s results for 


laminar flow. As Re* is increased the pumping action of large 


Figures 6, 7, and & show the variation of mixed-mean-to- 


Wall temperature differenes with Rayleigh number. Increasing 


tt 


Re* an order of magnitude approximately decreases @ m ai OPC ex 


of magnitude. Increasing Ra has no effect on Q. Increasing 
Pr has less effect on decreasing 0, than increasing Re*., 
Finally increasing F has the most effect on decreasing O rae 
As shown in Figure 8, Q,, passes through zero for F approxl- 
mately equal tol. If 0, is plotted against 1/F, as is shown 
Pa Eipure 95, tro aS Seen toate ivarie liom, in ha, Piscor Re> hes 
Mivule erlect On this Crossover point, ) Haliman’s curve for 
laminar flow, RazO is plotted for comparison 

When @y,=0, Nusco, according to equation D-3, Figure 10 
shows Nusselt number plotted against 1/F for Ra=0, Re*=103, anid 
Pe Or Ieard 10, 2.80 ini Alle discontinuity Occurs  wienever 
O a79- Krom equation D-3 it is “Seen that Nu=0 for Fel. . Hall= 
man's plot for laminar flow, Ra=0O is shown for comparison. 


- 
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Representative velocity and temperature Drei idtecry are Sno 
an Pigeres ll through el. 

The effect of varying Re* is shown in Figures 11 through 
16 for Pre10, Ra=625, and for various values of F. As Re* 
tnereases from 102 to 10", the velocity iy the center of ene 
tube becomes less positive and the velocity near the wall be- 
comes more positive. As Re* is increased, the temperature 
profiles flatten out and the slope becomes higher at the wall. 
Decreasing F, as shotm in Figure 21, flattens the velocity pro- 
file, whereas the temperature profile decreased considerably 
eiig  Dasses UNrOUgiaZerO auel=l, Vithe errecr Of 2nCrTeasiue 2h on 
the velocity and temperature profiles was not as is predicted 
Mima Pcmidim ! ., Ine aeulLen 15° Cue sLO the tincreascd values or 
Re* investigated in the turbulent flow problems. 

Figures 17 and 18 show the effect of various Prandtl nun- 
Ders 1 oT Re*=10°, Bootons and F=10. As Prandtl number increases, 
the velocity and temperature profiles decrease more in the cen- 
ter of the tube than near the wall, thus becoming more flat. 

The effect of a variation of Ra for Re*=10-, Pace 10 4 cane 
fa.5> 35 SAOun In Hieures. £9 and 20-> ss tha lis snercaccd trom 0 


ut 


BO 10 Stunereoule Very aipcle GCecreoascianetvc VeLlOci6., cane “cen 


perature Proliles at tne center of Che tube. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 


The effect of the parameters Re*, Pr, Ra, and F on Nu, C, 


ce, U,; and @ profiles is summarized in Table e. 


TABLE 2 


SUMMARY OF EFFECTS IN COMBINED 
FORCED AND FREE CONVECTION 


Changseéon Nu C 
for | Rat625 Rad 625 


@s 


increasing 


Re* increase increase increase 


Pie increase no effect decrease 


Ra no effect no effect increase 


30 


mn U 
Velocity @enperavure 
Dirolate Git erence 
Diol tac 


decrease flattens flattens 


Curve curve and 
makes 
approach 
Zero 

decrease lowers same as 
velocity Re* 

in cen- 

Cer OL 

tube, nak- 

ing it 

less pos- 

hieive 


Sliche increase decreases 
decrease in center in center 
ao RiLem Or Gone of tube 
Ra 


a, 2 





Table 2 (Continued) 


Change on Nu 


Meg 
Bo 
increasing 
Re* O 
ie 0 
Ra 0 
lca 
increasing 


Re* increase 


Pr increase 


ha TO eCflecy 


decreasing 
li slight 
increase 


Ra< 625 


increase 
no effect 


no effect 


increase 
no effect 


no effect 


no effect 


Ra>625 


increase 


deeneasec 


no effect 


increase 


increase 


Gecrease 


decrease 


{ad 
i 


Prn 


decrease 


aecrease 


no effect 


aeerease 


aecrease 


no effect 


increase 


U 
Profile 


flattens 
curve _ 


no effect 


no effect 


flattens 
curve 


flattens 
curve 


flattens 
curve 
slightly 


flattens 
curve 


p 


Profile 


negligible 
errect 


negligible 
effect 


no effect 


same as 
ag a 


same as 
teak 


gkayeelki vera iow Ms: 
Ciieser 


increase, 
passing 
through 
zero, Fal 





Conclusions 


ite Miimoes bh orOllLemon CieChien SUrTiICVentLly GlOose TOmpreylous 
known results to insure the validity of the program. 

ea The new expression used for the relative viscosity enabled 
tne universal velocity and temperature profiles to be approached 
almost exactly in pure forced-convéction turbulenv Went transfer 
problems with negligible volume heat sources. 

B. VOLUME eheatecourecs Had ner lie i bre teliceteon tnervclociTy 
profile, but had considerable effect on the temperature profile 
as would be expected from the development of the equations. 

4, The program was unable to solve problems involving high 
Rayleigh numbers and volume heat sources for low values of Re*. 
Dis increasing the various input parameters tend to flatten 

the velocity and temperature profiles as would be expected. 

or The effect of volume heat sources was reduced considerably 


py the dniiuence of turbulent Tlow. 
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Figure 2, Comparsion of Problem 429B and 397 Results 
with Universal Velocity and Temperature 
Profile 
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Figure 4, Variation of Nusselt Number 
with Rayleigh Number 
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PRESSURE DROP PARAMETER, C 


oe Eee te Bree es ae ee “et ee ee em ee eri te ee em POT COE LIE RL IE LE LE BG Et 


ik O LN I Ee A 0) etre nr ee 9 et ene: | e-em —s 








mash Snares 1D, el) BON Weel Cae 


1. 1.0 LO 10° 





RAYLEIGH NUMBER, Ra 


Figure 5. Variation of Pressure Drop Parameter 
With Rayleign Nunhepy 
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Figure 7. Variation of Negative Mean Tem- 
perature Difference with Rayleigh 
Number for F = 10 
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APPENDIX B 


THE EDDY DIERUSTVITY OF MOMENTUM 


An analytical study was made by Jackson! of the eddy 


Grrusivity Ol Momenvumn 1a vertical circular tube, wns 


2 ,. 4 Ge = Pee eee “, 
ratio of eddy diffusivity to kinematic viscosity, -Ef- , 
+’ 


Calied relative viscesity appeared in a reduced form of cvhe 


Navier-Stokes equation. An empirical fit of experimental data 


was employed to determine the dimensionless velocity gradient 


appearing in the expression for fe eel 


VS a8 


tout 


The result of His Study used in the present analys 
Pane LOLMLOWING SeOlUL CXR Oress 0OmS {Or Veariatien Cf relacive 


viscosity with dimensionless radius j 


m 
= 
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The results commare favorably with the experimental and 
analytical results of others and annuears to be an improvement 
over the £m, GAO Seale eee ap cece: = budd ca Teese. Tee 


sults appear as equations (29), (30), and (31) in the present 


analysi 





APPENDIX C 


THE LONGITUDINAL TEMPERATURE CRADLENT 


The heat balance for the steady-state condition in our 
problem states tnat the rate at which heat is added to the 
fluid is equal to the rate at which heat is etved down-~ 
stream in the fluid. 


Tn symbolic form (see Figure C-1), 


“eddy oq phd - fe fuemrdec dtd; (C-1) 
} dy bw . [Gn f a dx ‘ 


where q,' is taken positive if heat is being 


& removed, 
for ty give 

at ba a. | (c-2) 
Ox CP Cry, Ve furdr 


The mean velocity is 


solving 
equecion (C-1) 


defined as follows: 


D,. Oto 
Un, a fue U2irrar | _ as LA rar (c=2) 
— 19 D ee op On ie feo ras > 
from which 
LY; 
Z furdr — Uys (c~h) 
0 § 


Equation (C-2) can therefore be written as 


Oe pe 2 a. 


VAs 


we (C~5) 
3x P rr Cio yn [) 


Since all quantities on the right~-han 


side of equation (C-5) 
are constant in our problem, 





oy = EH (a constant) (C-6) 





FIGURE C-1 


Heat Balance for an Element of Fluid 


Therefore, the temperature gradient in the direction of 


flow is a constant independent of radial location. 


61 





APPENDIX D 


NUSSELT MUMBER 


Nusselt riumber Nu is a useful expression in comparing 
our problem results. it is Function ef the tempcrature grad- 
menG ay the Wall (CO a constant average Genperature eradient?, 


Considering anes Velaree Oi a walt Lercth of tubes 


[= rb a q PD = Cras CaF] zp (D-1) 


where gq," is taken positive if heat is being removed? , Sub- 
stituting in equation (D-1) the equations 


Cy = h Ob 
LEK Brn 
yD 
; be AD 
SLves a ne =e GDe<2) 
Ge Um Co Ff K ae oem | 


Substituting the definition of F in equation ({D-2) and rear- 


and 


2. 
on 
! 





ranging terms gives 


Nu = hD/k = = (1-F) (D-3) 
Pin 


From this equation, as Oy, appre Baca. 2ZCrone eo sae 


proach infinity; and when F = 1, Nusselt number will equal zero. 
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APPENDIX Is 


REYNOLDS NUMBER 


Reynold number Re is the ratio of the inertia forces to 


~ 


frictional forces ari is e measure of turbuienc 


Considering a dyrnamie force balance on an element of 


fluid, the following expression can be developed: 


we ede - 5 4 De J (E-1) 


ox dx ep Sis 
where { is for upward flow and - is for downward flow, 
From equation (13) the pressure-drop parameter C is 


C= “Ore (GE + v Gan (E-2) 


32M Un 


and substituting equation (E-1) in equation (E-2) gives 


(tal Sees (B-3) 
ae SEM Uy 
Using eee. 
Fw = (u") eae 
Cw - km = Ow BOm | 
Bin = ee DO (eae 


and the parca of Re, Re*, Ra, and F, we obtain 


nat 


s t ta ee Ra va 
= Oe oo (B-) 


Solving for Re sae 


eae (3-5) 


CR a = ee 


which is equation (35) in the Analysis. 
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APPENDIX F 


VOLUME HEAT SOURCES 


The term a''™ accounts for the volume heat sources in 


our problems, Since only uniform volume neat sources are 
treated, g''!! is independent of direction, 


_— 


The nondimensional volume heat source parameter is 


= C LA wn /e _— 
Bess ies ee (F 1) 


as defined by Hallman. The physical significance of F can 

be GCeseribedias Une ratio Of the tnermal enersy convected 
downstream, per unit volume to the heat generated in the fluid 
per unit volume. Considering a heat balance on a unit length 


of pipe with constant wall heat addition 


i] 


TT e | . ; 
] 4 ~ Jw atyineeta te. (F-2) 


and introducing equation (F-1) gives 


i] 
ae ee = ban jae : (F-3) 


The meaning of the various values of F can be seen from 
equation (F-3). When qi = 0, which corresponds to an insulated 
Lupe Wall then Gos i pneiWwalls should be notver than the 


bulk of the fluid because heat is conveeted away from this 
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regi@n less rapidly than in the center, and because of uniform 
volume heat source peneration. Fluid velecities should there- 
fore tend to be higher in the vicinity of the wall than near 
mnie Cenver line, 

When F< 1, it means that heat is being removed at the 
walls because more heat is being generated internally than is 
convected downstream, The velocity at the wall should be slow- 
er than the velocity at the center line in this case since heat 
wo  Oeine removed and tne density of the fluid as" hieher av. tie 
wall. 

When F>1, it means that heat is being added at the walls 
because heat is being convected downstream faster than it is 
being generated. The velocity in this case should rise faster 
near the walls than near the center line since it is hotter at 
the walls. 

R= Owe occur melencver there (a no net ehrouci=2 lo 
(un = 0), no flow (u = 0), or constant wall temperature (Azx0). 
An investigation is being conducted on the general problem 
WiEh a constant wall temoerature by Pagel’, 

FS eCo) Would corres pond. UG ne volume heat sources.” An 
investigation is being conducted on this general problem by 


Grirfinl®, 
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